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We calculate the magnetic-field and temperature dependence of all quantum corrections to the 
ensemble-averaged conductance of a network of quantum dots. We consider the limit that the 
dimensionless conductance of the network is large, so that the quantum corrections are small in 
comparison to the leading, classical contribution to the conductance. For a quantum dot network 
the conductance and its quantum corrections can be expressed solely in terms of the conductances 
and form factors of the contacts and the capacitances of the quantum dots. In particular, we calculate 
the temperature dependence of the weak localization correction and show that it is described by an 
effective dephasing rate proportional to temperature. 
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I. INTRODUCTION 

The low temperature conductivity of disordered met- 
als or semiconductors is dominated by the elastic scat- 
tering of electrons off impurities and defects. While the 
conductivity is determined by Drudc-Boltzmann theory 
for not too low temperatures, quantum corrections to 
the conductivity become important at temperatures low 
enough that the electronic phase remains well defined 
over distances large in comparison to the elastic mean 
free path.ErHa One usually distinguishes two quantum 
corrections, the we&kJacalization correction and the in- 
teraction correctionJatra The former is caused by the con- 
structive interference of electrons traveling along time- 
reversed paths, whereas the interaction correction can be 
understood in terms of resonant scattering off Friedel os- 
cillations near impurities .oH 

Although they are small in comparison to the Drude 
conductivity, the quantum corrections are important be- 
cause they strongly depend on temperature and an ap- 
plied magnetic field, whereas the Drude conductivity 
does not (as long as impurity scattering is the domi- 
nant source of scattering). Theoretically, the tempera- 
ture and magnetic-field dependences of the corrections 
can be expressed in terms of the sample's diffusion con- 
stant (or, equivalently, the elastic mean free path) , which 
can be obtained independently from a measurement of 
the Drude conductivity. The availability of quantitative 
theoretical predictions has made a detailed., comparison 
between theory and experiment possible. 

The same quantum corrections also exist for a 'quan- 
tum dot', a conductor coupled to electron reservoirs via 
artificial constrictions {e.g., tunnel barriers or point con- 
tacts), such that the conductance of the device is domi- 
nated by the contacts and not by scattering off impurities 
or defects inside the sample. The latter condition is sat- 
isfied if the product E^v of the dot's 'Thouless energy' 
and its density of states is much larger than the dimen- 
sionless conductance of the contacts connecting the dot 



to source and drain reservoirs. (The Thouless energy is 
the inverse of the time needed for ergodic exploration of 
the quantum dot.) 

In this article we consider 'open' quantum dots, which 
have contact conductances larger than the conductance 
quantum e 2 /h. Because transport through a quantum 
dot is dominated by the contacts, it is described by the 
sample's conductance, not its conductivity. The quan- 
tum corrections then pertain to the conductance after 
averaging over an ensemble of quantum dots that differ, 
e.g., in their shape or precise impurity configuration. 

While the magnetic-field dependence of quantum cor- 
rections to the ensemble averaged-conductance is in ap- 
parent agreement with the theory,E3 the situation regard- 
ing the temperature dependence is more complicated and 
no good agreement has been reported to date. Theoret- 
ically, the temperature dependence of the weak localiza- 
tion correction to the conductance of a quantum dot is 
described by means of a 'dephasing rate' 7^. For a quan- 
tum dot, one expects 



70 = cT 2 /E. 



Th 1 



(1) 



where T is the temperature and c is a numerical constant 
that depends on the dot's size and shape. E3ll3't£l The pro- 
portionality constant c can not be measured indepen- 
dently, however, which is an important difference with 
the case of a diffusive conductor. The absence of a sepa- 
rate method to determine this constant poses a significant 
difficulty when comparing theory and experiment. A sec- 
ond difficulty is the lack of a direct theory of the tempera- 
ture dependence of weak localization. Instead, the avail- 



xiiitions employ a phenomenological 
1 1 ~ 5 ' and match the dephasing rate 



able theoreiij. 
descriptionEM 

to Eq. (Q), from which the temperature dependence of 
weak localization can be obtained. 

In this article, we study the quantum corrections to the 
conductance in a network of quantum dots or "quantum 
circuit" .E3 (See Fig. |l| for an example of a quantum dot 
network with A/d = 3 dots.) Replacing a single quan- 



2 




FIG. 1: An example of a quantum dot network with Ad = 3 
quantum dots. The conductance of the network is dominated 
by the conductances of the contacts between the dots. We 
assume that all dots in the network are 'open', i.e., all contact 
conductances are much larger than the conductance quantum 
e 2 /h. 



turn dot by a network solves both difficulties mentioned 
above: A quantum dot network allows a calculation of 
the complete temperature dependence of the quantum 
corrections to the conductance without the need of an in- 
termediate step involving a phcnomcnological dephasing 
rate and without parameters that can not be measured 
independently. The relevant parameters in a quantum 
dot network are the conductances and form factors of 
the contacts in— the network and the capacitances of the 
quantum dotsEH 

Our main result is an expression for the ensemble av- 
erage of the dimensionless conductance 



n d s e 2 



(2) 



where d s — 1 or 2 in the absence or presence of spin 
degeneracy, respectively. The result becomes exact in 
the limit that the contact conductances are much larger 
than the conductance quantum e 2 /h, 



(g)=9 cl + Sg WL + Sg m ^ + 



(3) 



Here g cl is the 'classical' conductance one obtains from 
Drude-Boltzmann theory, while 5g WL , Sg™ 1 ' 1 , and 5g mt ' 2 
are three quantum corrections to (g). Explicit expres- 
sions for g cl and the three quantum corrections in terms 
of the contact conductances and the capacitances of the 
quantum dots in the network, as well as the precise con- 
ditions for the validity of Eq. ([|) will be given in Sec. Q 
below. The correction 5g WL is the weak localization cor- 
rection. It is the only quantum correction that is affected 
by the application of a magnetic field. The remaining two 
corrections arise from electron-electron interactions. The 
first interaction correction Sg 111 *' 1 represents a non-local 
correction to the conductance t hat earis ts for networks of 
two or more quantum dots onlyr 9 IP°H 31 l It is the counter- 
part of the Altshuler-Aronov correction in the theory of 



disordered conductors. The second correction, 5g lnt ' 2 , de- 
scribes the renormalization of the contact conductances 
by the interactions. It is usually referred to as (dynami- 
cal) Coulomb blockade, an effect that is well-known from 
the theory of transport through tunnel junctions in se- 
ries with a high^mpedancfi or quantum .dots, with tun- 
neling contactsMHmmMBBm.By Its coun- 
terpart in the theory of disordered conductors is the 
Altshuler-Aronov correction to the tunneling density of 
states.Ej 

The fact that the temperature dependence of quantum 
corrections in a quantum dot network does not depend 
on details of individual dots has its origin in the different 
form of the relevant electron-electron interaction modes 
in a quantum dot network and in a single dot. In a single 
quantum dot, the dominant contribution to the electron- 
electron interaction is the uniform mode, the strength of 
which is set by the dot's capacitance. Apart from a possi- 
ble renormalization of the contact conductances, 5g lnt ' 2 , 
the uniform mode has no effect., n.n the quantum correc- 
tion to the dot's conductance.OC3a'C3 In particular, the 
weak localization correction Sg WL is unaffected by the in- 
teraction and the non-local interaction correction 8g mt l 
vanishes. Instead, electron-electron interactions deter- 
mine (>g WL and Sg lnt ' 1 in a single quantum dot through 
sub-dominant non-uniform interaction modes, rWjlaich are 
known to depend on the precise sample details.tZlEZI For a 
quantum dot network, on the other hand, there exist in- 
teraction modes that are uniform inside each dot but not 
across the full network. With such interaction modes, all 
three interaction corrections Sg wh , Sg" 1 *' 1 , and 8g mt ' 2 are 
generically nonzero and temperature dependent. More- 
over, because these modes are uniform inside each quan- 
tum dot, their properties depend on the contacts between 
the dots and on the dot capacitances only, not on the 
precise geometry of each dot separately. It is this essen- 
tial feature that makes a quantum dot network an ideal 
paradigm for studying the effect of electron-electron in- 
teractions on quantum transport in finite-size systems. 

Separate aspects of the problem we address here have 
been considered before. Weak localization in single quan- 
tum dots withaut. interactions has been studied by var- 
ious authorspH-BEffia-BBH as well as the ef- 
feet of the uniform interaction mode on the conduc- 
tances of the cput s .ctH-aMM^a&-i^^Qi.^ the elec- 
tron reservoirs.BByee-ElBEMjy>BBy (See Ref. 
^ for a discussion of a comparable effect involving 
spin-dependent interactions in the quantum dot.) Also, 
while it is known that the uniform interaction mode 
has no effect on weak localization because a spatially 
uniform fluctuating potential affects . pha ses of time- 
reversed trajectories in the same way,c3c3 the uniform 
interaction mode can suppress interference contributions 
to other obser*ables-i if the quantum dot is part of an 
interferometer B'Ej'Ey 

Weak localization in networks of quantum dots, but 
without interactions, was considered, by Argaman for 
dots connected by ideal contacts,EirE3 and by Cam- 
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pagnanOj-apd Nazarov for dots connected by arbitrary 
contacts.EjGolubev and Zaikin calculated the interaction 
corrections-j^g 111 *' 1 and <5g lnt ' 2 for a linear array of quan- 
tum dotsj23 as well as the weak localization correction 
for non-interacting electrons (but with a phcnomcnolog- 
ical dephasing rate).E3 In a recent publication, the same 
authors also considered the full temperature dependence 
of weak localization in the special case of a double quan- 
tum dot (a network, with Ad = 2 quantum dots) with 
tunneling contacts j£j and reported that electron-electron 
interactions suppress weak localization even at zero tem- 
perature, a conclusion that contradicts the common wis- 
dom that there is no dcphasing n ftom electron-electron 
interactions at zero temperature. ErEI 



Weak localization and interaction corrections have 
also been, considered for networks of diffusive metallic 
wires.oE3 Large arrays of quantum dots connected by 
tunneling contacts further appear in the study of granu- 
lar metals.E3 Beloborodov and coworkers considered the 
interactipn-^QXijcctions 6g lntl and 8g mt2 for a granular 
metalJ^OQOEJ but accounted for weak localization 
and its temperature dependence only via a phcnomeno- 
logical dephasing rate and a renormalizcd diffusion con- 
stant. A microscopic theory of the temperature depen- 
dence of weak localization in granular metals was-given 
by Blanter et al. in the high temperature limit £il Our 
present analysis (as well as that of Ref. [30]) is for con- 
tacts of arbitrary transparency and contains contribu- 
tions to weak localization and to the interaction correc- 
tion to the conductance that are absent in a network 
where all contacts are tunneling contacts. Our results 
agree with the literature wherever applicable, except for 
the zero-temperature limit of the weak localization cor- 
rection Sg WL , where we find that weak localization is 
unaffected by electron-electron interactions, in contrast 
to Ref. Ejl 



The remainder of our article is organized as follows. In 
Sec. [jl] we introduce the relevant parameters needed to 
describe the quantum dot network, formulate our main 
assumptions, and present our main result, an expression 
for the ensemble-averaged conductance and its quantum 



corrections. In Sec. Ill we motivate our result for the 



temperature dependence of the weak localization correc- 
tion using semiclassical arguments. In Sec. IV we then 



turn to a fully quantum mechanical calculation of the 
conductance and its quantum corrections using random 
matrix theory. We specialize to the simplest network, a 
double quantum dot, in Sec. [v| and discuss the origin of 
the difference between our result and Ref. ^33] for the zero- 
temperature limit of weak localization. We conclude in 



II. DEFINITION OF THE PROBLEM AND 
MAIN RESULTS 

A. Network of quantum dots 

We consider a network of Ad quantum dots, coupled 
to two electron reservoirs. A schematic drawing of a net- 
work is shown in Fig. In this section we introduce the 
relevant parameters to describe the quantum dot network 
and summarize our main results. 

The quantum dots are connected to each other and to 
source and drain electron reservoirs via point contacts. 
The dots will be labeled by an index i = 1, . . . , Ad; the 
reservoirs are labeled by the index a = 1,2. The contact 
between dots i and j is described by its dimensionless 
conductance (per spin direction) and its form factor 
fij. Both gij and fij are defined in terms of the trans- 
mission matrix fy of the contact, 



>r,j 



fij — tr (tijt^j) 



t ^2 



(4) 



Form factors are related to Fano factors /? often encoun- 



tered in the literature via /3j. 



(ffij - fij)l9i. 



The 



dimensionless conductances and form factors are sym- 
metric, g^ = g jt and f v] = f jl7 ij = l,...,Ab- Spin 
degeneracy will be explicitly taken into account via the 
parameter d s = 1,2. 

Similarly, the contacts between the ith quantum dot 
and reservoir a, a = 1,2, are described by a dimension- 
less conductance g' ia = g' ai and a form factor f[ a = f' ai , 
which are related to the transmission matrix t' ia of these 
contacts as 



tri' i' f 



f - tr (t' t n ) 2 

Jia — L1 \ L ia L, ia) 



(5) 



For ballistic contacts one has / = g; for tunneling con- 
tacts one has / « j. Throughout we assume that all 
conductances are large, 



9ij, #L:> !j h3 = !,■•• ,A/b, 



1,2. 



(G) 



(One may replace this condition by the less strict re- 
quirement that each quantum dot be well connected to 
one of the two reservoirs, such that the regime of strong 
Coulomb blockade is avoided.) For future use, we arrange 
the conductances and form factors in A/"d x Ad matrices 
g and / with elements 



9ij = 
fij = 



a=l 9aj 

~9ij 

V 2 f ■ 

~fij 



Z^fe^i 9ik 1 — J, 

i + 3, 

12k^i f ik * = 3l 

i ^ 3- 



(7) 
(8) 



Sec. VI 



The quantum dots are assumed to be disordered or 
ballistic-chaotic, with density of states v i per spin degree 
of freedom and Thouless energy i?Th,i, i = 1,...,Ad- 
The Thouless energy -Btii.i = ^/lerg.i, where T ergi i is the 
time for ergodic exploration of the ith quantum dot. If 
the electron motion is diffusive inside each quantum dot 
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with diffusion constant D, £/rh.i ~ D/L 2 where Li is the 
linear size of dot i. (Our definition, while common in the 
literature, differs from some references where £?Th,-t is the 
inverse of the dot's dwell time.) We assume 

E-r^iVi > gu, i = 1, . . . ,7V D , (9) 

so that random matrix theory can be used to describe 
the electronic states in the quantum dot network. An 
external magnetic field is described by means of the di- 
mensionless numbers 

gu,i = E Thti Vi-^, i = l, ...,Nt>, (10) 

where $i is the magnetic flux through the ith quantum 
dot and $o — hc/e is the flux quantum. In order to 
simplify the notation, we arrange the densities of states 
Vi and the parameters <?H,i m diagonal A/o-dimensional 
matrices v and gu, 

i>ij = i>i5ij, (gn)ij = gn,Aj, i,j = 1, • • • ,-A/d- (11) 

Corrections to the conductance that depend on the mag- 
netic field will only be relevant where gu,i is of order gu 
or less, otherwise they will be fully suppressed. In that 
parameter range, the flux through the insulating regions 
between the quantum dots is much smaller than $o, so 
that the corresponding Aharonov-Bohm phases can be 
neglected. 

The inequality (^) also implies that the ,-eJectron- 
electron interaction in each dot is well screened.eH Hence, 
the electron-electron interaction couples to the total 
charge qt = erij of each dot only. Such an interaction 
is described by means of capacitances for the capac- 
itive coupling between dots (if i j) and for each dot's 
self-capacitance (if i = j). Again, we arrange the capac- 
itances into an Ao-dimensional matrix C , 

6ij = [Tt^ ik i^j, (12) 

For metallic dots, one has the inequality 

Cn/e^^Ui, i = l,...,Nh- (13) 



B. Quantum corrections to the conductance 

Our main result is a calculation of the ensemble- 
averaged conductance (G) = (d s e 2 /h)(g) of the quantum 
dot network as a function of temperature, 

(g)=g cl + Sg WL + 6g^ 1 +Sg^ 2 , 

where g cl is the classical conductance of the network and 
<5<? WL , Sg 11 ^' 1 , and 5g mt ' 2 are corrections. The average 
conductance is calculated using the following limiting 
procedure for the parameters of the network: 



1. We first take the limit (0) needed for the applica- 
bility of random matrix theory, while keeping the 
ratios Vi/vj and T/z^, as well as the <?h,i fixed, 
i,j = 1, . . . ,Af B . 

2. We then take the limit ([|) of large contact conduc- 
tances, while keeping the ratios gij/gik, gij/gn,i, 
and g i3 /g' ia fixed, i, j, k = 1, . . . ,Ad, a = 1,2. 

3. Finally, we simplify our results using the inequality 
(p^|), if possible. 

In all three limiting steps, the number A/d of dots in the 
network is kept constant. Keeping the ratio T jvi fixed in 
the first limiting step eliminates interaction corrections 
from non-uniform interaction modes inside the quantum 
dots, see Eq. ([j]) above. In the second limiting step, we do 
not make any assumptions about the temperature, thus 
allowing for the full range of temperature-dependent ef- 
fects that can be described within random matrix theory. 
We note that, while the classical conductance g cl diverges 
in this limiting procedure, this divergence does not affect 
the temperature or magnetic-field dependence of (g) be- 
cause g cl does not depend on temperature or magnetic 
field. Corrections not included in Eq. (|J) are either small 
in the limit (ph of large contact conductances or small 
in the limit (0) used to justify the use of random matrix 
theory. 

The leading term g cl in Eq. (||) reads 

A/d 

= si.fl.TVa, (14) 

where, in the second line of Eq. ([l4|), we have written "•" 
to denote indices in adjacent factors that are summed 
over as in matrix multiplication. [Compare with the first 
line of Eq. (p)).] This shorthand notation will be em- 
ployed throughout the text. 

The correction Sg WL is the weak localization correc- 
tion to the ensemble-averaged conductance. It can be 
distinguished from the remaining two corrections Jg 111 *' 1 
and Sg mt ' 2 because 5g WL depends on an applied magnetic 
field whereas Sgf* and 5g lnt > 2 do not. We find 

6 9 WL = ->Y. Ziia'iXg^-r^ig-hfgjh'.i 

A/d 

i=i 

A/d 

+ E £ " g'l-^ifi-g' 1 ^ - fa) 

i=l 

Ad 

- E /« % Vt V 2 . (1 5 ) 
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where the Ad x Ad matrix c is the counterpart of the 
"Cooperon" in the theory of weak localization in disor- 
dered conductors. For the quantum dot network, c reads 

^ = E^r( r + r H + r0)r fc :, feI (is) 
fc=i 



where T, Th, and are rank- four tensors, 



2nhu 



1 - i 1 r ~ 

ffifcdj; + — - — dikgji 
Zirtii/j 



2~K%Vj 



(F<l>)ik,jt = -JJ-idu 1 +.9j/ - ^g^SikSji. (17) 

The terms Th and describe the suppression of weak 
localization by a magnetic field and electron-electron in- 
teractions, respectively. In the limit of low temperatures 



= and Eq. (16) simplifies to 

hj = {g + m)T/- 



(18) 



For high temperatures (T^ajj diverges [other elements 
are zero because of the Kronecker deltas in Eq. (\n\)], 
except for the diagonal elements with i — j. Hence, one 
finds 



c I7 - = {g + gn)a , 



(19) 



where gfj is the diagonal part of the matrix g, gf- = (jij Sij . 
This is the contribution to the weak localization cor- 
rection that arises from self-returning electron trajec- 
tories that reside inside one quantum dot only and, 
hence, are unaffected by dephasing from electron-electron 
interactions .Lil 

The first interaction correction <5g mt,:L is 



27r f ( d io \ ^ D 

' a, 13=1 k,l=l 



Vfj 



d B J \duj 

x(g- 2nii>u))~l(g - 2nii>uj)ki(g - 2niDuj)^lg 1 .(g7 a 1 - gT^)^ 1 - ~g~ X )g'.i 



(20) 



The second interaction correction <5<? mt ' 2 represents the 
renormalization of the conductances between the quan- 
tum dots and between the dots and the reservoirs as a 
result of the electron-electron interactions, 



Sg 



Ar> 2 „ Ad n 

int,2 d 9d sJ , <J9c\ 



The interaction corrections 8g' ia and Sgij exist for non- 
ideal contacts with f io < g i} , f' ia < g' ia only, i,j = 
l,...,Ab, a = 1,2, 



&g'aj = -(g'aj - fad 



$9jk = ~(9jk - fjk) 







coth 














2T 


du> 




coth 


LU 














2T 



x Re (Szjj + 5zkk - 2Szjk), 



(22) 



where 5z is the difference of the network's dimensionlcss 
impedance matrices with and without interactions, 

Sz = (d s g — 2itiujC / e 2 ) 1 — (d s g — 2iriu}d s i')~ 1 . (23) 

The interaction correction Jg 111 *' 1 was obtained previ- 
ously by pplubev and Zaikin for a linear array of quan- 
tum dots£j and by Beloborodov et al. in the context of a 



granular metal.El It is the counterpart of the Altshuler- 
Aronov correction in disordered metals, where it arises 
from the diffusive dynamics of the electrons. Although 
the electron dynamics is not diffusive in a quantum dot 
network, it is non-ergodic, which is sufficient for this in- 
teraction correction to appear. (The exception is a quan- 
tum dot network consisting of a single quantum dot only, 
for which the electron motion is ergodic. Indeed, one ver- 
ifies that iSo 1 " 4,1 = if Ad = 1, m agreement with Refs. 
|30|,^1 43 5|.) A semiclassical calculation of Sg™ 1 ' 1 for the 
special case of a double quantum dot with ballistic con- 
tacts can be found in Ref. |3l[ 

For the case of a single quantum dot, the renormal- 
ization of the contact conductances Sg lnt ' 2 or "dynami- 
cal Coulomb blockade" was obtained previously in Refs. 
||,|9|j40[^l|,||,|3[g. The renormalization of the contact 
conductances in the quantum dot network is essentially 
the same as in the case of a single quantum dot or a single 
tunnel junction coupled to a high-impedance electrical 
environment — in both cases the change of the contact 
conductance is proportional to the factor (g — /) — , the 
only difference being that the impedance z is replaced 
by the impedance matrix z in the case of the quantum 
dot networkEj The same conclusion was reached for the 
interaction correction in an array of quantum dots with 
tunneling contacts ip .thft-context of transport through a 
granular metalMfflM 
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Equations (||)-(|2^) provide a general solution for the 
ensemble-averaged conductance and its quantum correc- 
tions in an arbitrary quantum dot network for arbitrary 
temperature. These expressions can be simplified only 
by specializing to a particular quantum dot network. In 
Sec. [V| we analyze these expressions for the case of a dou- 
ble quantum dot, a network consisting of two quantum 
dots. 

Although it is not possible to proceed quantitatively 
without specializing to a particular network, we can com- 
pare the sizes of these three quantum corrections and 
their typical temperature dependences. For the limiting 
procedure taken here — see the discussion following Eq. 
(||) — , the relevant temperature iSCale for dephasing of 
the weak localization correction isO 



where 



T,/, = hmax(g,gn)/TD, 



td ~ hv/g 



(24) 



(25) 



is the typical dwell time for the network. (Here g and 
gn are shorthand notations for typical values of gij or 
gu,i in the network, respectively.) For the interaction 
corrections Sg lattl and 5g lnt,2 j the relevant temperature 
scales are Ti/td and the inverse charge relaxation time 



h/r c ~ e 2 g/hC. 



(26) 



(In a more precise analysis one needs to identify Ad dwell 
times and A/d charge relaxation times for a network con- 
sisting of Ad quantum dots, see Sec. |v| for an explicit 
calculation for Ad = 2.) Since, typically, C/e 2 <Si v, 
the charge relaxation time and the dwell time satisfy the 
inequality 



(27) 



With these definitions, we find the order of magnitude 
of the weak localization correction Sg wh to be 



Sg 



WL 



s 9wl 



max(l, T/T^) 



(28) 



where Sg^ L an d ^3wl are constants of order 
min(l, g/gn)- Similarly, for interaction corrections we 
find 



Sg 
Sg" 



int.l 



min(l, %/Ttq), 

ln[max(r c T/?i, r c /r D ) 
U/Tt c 



if T«n/r c , 
if T>ft/r c , 



(29) 
(30) 



independent of the magnetic field. All three quantum 
corrections need to be taken into account for a complete 
description of the temperature and magnetic-field depen- 
dence of the conductance of a quantum dot network. In 
particular, in order to correctly describe the temperature 
dependence of (g) for T < %/rr>, Sg 1 ^' 1 can not be ne- 
glected with respect to 8g lnt ' 2 , in spite of the fact that 




FIG. 2: Schematic drawing of a trajectory a and its time- 
reversed a that contribute to the Cooperon propagator c. 



5g mt ' 2 is larger than Sg 1 ^' 1 by (at least) a large logarith- 
mic factor 1ii(td/t c ). 

The temperature dependence (|8|) implies a dephasing 
rate that is linear in temperature. A linear tempera- 
ture dependence of the dephasing rate was obtained pre- 
viously-by Blanter et al. in the context of a granular 
metaljlil and by Seelig and Buttiker for a single Quan- 
tum dot embedded in one arm of an interferometer .E3 In 
both cases, the linear temperature dependence of the de- 
phasing rate arose because the fluctuations of the electric 
potential can be considered classical, similar to the situa- 
tion encountered in one-dimensional and two-dimensional 
disordered conductors £3 As we discuss in the following 
sections, the same mechanism is responsible for the lin- 
ear temperature dependence of the dephasing rate in the 
quantum dot network. 

In the next section we describe a semiclassical deriva- 
tion of the weak localization correction and its temper- 
ature dependence, Eq. ( |l5| ) above. A full quantum me- 
chanical calculation of all three corrections to the con- 



ductance is given in Sec. |V|. We apply the general results 
presented here to the specific case of a double quantum 
dot in Sec. [v|. 



III. WEAK LOCALIZATION: SEMICLASSICAL 
CONSIDERATIONS 

In this section, we give a semiclassical argument for 
the temperature dependence of the weak localization cor- 
rection to the conductance of a quantum dot network. 
These arguments provide a semiclassical interpretation 
of the fully quantum mechanical calculations of the next 
section. 

Weak localization appears because of constructive in- 
terference of time-reversed trajectories. This interference 
leads to a small increase to the probability P re t that an 
electron returns rto,its point of origin. Following the stan- 
dard arguments p'H P ret is calculated as a square of the re- 
turn amplitude which, in turn, is written as a sum of am- 
plitudes A a over all returning paths a. (These paths are 
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classical paths in ballistic conductors,E3K3 and quantum 
diffractivc paths in conductors with impurity scattering.) 
The quantum correction to P rot then follows from inter- 
ference between a path a and its time-reversed a. Since 
the length of the self-returning path is arbitrary, the weak 
localization correction to the dc conductance is propor- 
tional to the time integral of the interference correction to 
the return probability known as the "Coocjeron" in the 
diagrammatic theory of weak localizationBS The coun- 
terpart of the Cooperon for the quantum dot network is 
the quantity 



1 



(2'nfi) 2 ViV. j 



^ ' A a (A a ) 



(31) 



where the sum is over all trajectories a that originate in 
dot j and end in dot i and a is the time-reversed of a, 
see Fig. ||. [Note that the return probability involves the 
diagonal elements Cu of the Cooperon matrix only. We 
have included non-diagonal elements in Eq. (^lj) above in 
view of the discussion of interaction effects below. Non- 
diagonal elements with i and j in adjacent dots also 
appear for the description of weak localization in a net- 
work of quantum dots with tunneling contacts, see Eq. 
( |l5| ) above.] 

^At zero temperature and without a magnetic field, 
A a = A a . We may then calculate using that \A a \ 2 is 
the probability that an electron propagates along trajec- 
tory a. Hence 



1 r°° 



(32) 



where Pij(r) is the probability that an electron in dot j 
is found in dot i after time r. In Eq. ( [32] ) we canceled a 
factor 2-KTiVj in the denominator against the phase space 
volume of the jth quantum dot. For a quantum dot net- 
work, Pij(r) can be expressed in terms of a rate matrix 

7: 



P ij (r) = (e-^) ij) j = g/(2whu). 
Integrating over time, we then find 



(33) 



(34) 



The interference between a path a and its time- 
reversed is suppressed if time-reversal symmetry is bro- 
ken by a magnetic field, because a magnetic field changes 
the phases of A a and A s in opposite ways. Interference is 
also suppressed because of electron-electron interactions 
at a finite temperature. Interactions cause the electrons 
to experience a time-dependent potential (f>(r , t) , which 
modifies the phase of A a and A a in different ways if the 
trajectories a and a are in different dots at the same 
time t.i3 For a network of quantum dots, the fluctuating 
potential <j) is uniform inside each dot, so that we can 
write 4>{j,t), where j — 1, . . . , A/d is the index labeling 
the quantum date, in the network. For each amplitude 
A n one then hasE3 




dr 



+ 



t < 2dr 



dr 



t - 2dr 



FIG. 3: Calculation of the Cooperon propagator for a net- 
work of quantum dots. A trajectory a originating in dot j and 
ending in dot i and duration t is separated into two segments 
of duration dr and a remaining segment of duration t — 2dr 
if 2dr < t. A self-consistent equation for Cij is obtained by 
considering the combined effect of escape, the magnetic field, 
and the fluctuating potential to first order in dr. 



where t a is the duration of the path a, j a (t) the index of 
the quantum dot corresponding to the position of path a 
at time t, and A a [0] the return amplitude in the absence 
of the potential <f>. 

For a quantum dot network, one may consider as a 
classical fluctuating potential. (This will be verified in 
the exact quantum mechanical calculation of Sec. IV I 
below.) Its fluctuations are given by the fluctuation 
dissipation relationjlj 



(36) 

where the response function Lfj (cu) describes the (linear) 
change S(f>i/e of the electric potential in the ith quantum 
dot to a change Sqj = eSnj of the charge in the jth. 
quantum dot, 



84> t (lo) = —L^(u))5rij(u)). 
For the quantum dot network, one has 

LfAu) = - \(j/e 2 + d^V- 1 - 2Tiiujg- 1 )- 1 



(37) 



(38) 



where the matrices C, v, and g were defined in Sec. || 
above. Typically, Cu/e 2 Vi, gu/\w\, and we can re- 
place Eq. (11) by 



(39) 



(35) 



Using this expression for Lf-(uj), we find that Eq. ( |3^ ) 
simplifies to 

m,t)d>{j,t')) = ^L~g-U{t-t'). (40) 



In order to find the effect of the fluctuating potential 
on the Cooperon propagator cy, we separate the contri- 
butions from trajectories a of duration t a smaller and 
larger than 2g?t, where dr is a time interval sufficiently 
short that the net phase shift from the fluctuating po- 
tential in the exponent in Eq. ( |35| ) is small, see Fig. ||. 
We also take dr much shorter than the dwell time in a 
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single quantum dot, so that Py(dr) 



jijdr, see 



Eq. ( |33| ) above. For trajectories of duration t a > 2dr 
we consider the initial and final segment of duration dr 
separately. Recognizing that the contribution from the 
intermediate segments of duration t a — 2dr can again be 
expressed in terms of c, and using Eq. (|4^) to average 
over the fluctuating potentials, we then find 

2dr „ 



+ (5ifc - lkidT)(5ji - ^ijdr)c k i 



k,l=l 
A'd 



k,l=l 



-J — "ij 

Trnisj 



(r- 



^4>)ik, 3 iCkidT, (41) 



up to corrections of order dr 2 . Here 



7H,ij 

74w 



4ttT 



dJi 



(fl 



(42) 



andr,fej; = 7w5ji+6'ijfe7y ) (rn)i^ij = 7H,ifc^7+^fe7H,j/, 
(r^Jifc.j'i = %,ijSikSji, cf. Eq. (O) above. Solving this 
equation for c, we arrive at Eq. ( |16[ ) of the previous sec- 
tion. 

It is worth while to point out that the temperature de- 
pendence of weak localization is caused by processes that 
involve the exchange of energy quanta small in compar- 
ison to the temperature. Such processes are commonly 
referred to as "dephasing" , in contrast to more general 
inelastic processes which lead to a broadening of the elec- 
tronic distribution functionoEl In this sense, interaction 
effects in the quantum dot network differ from those in 
a single quantum dot, where weak localization is sup- 
pressed bjj inelastic processes that involve a large energy 
transfer O'ta Indeed, the characteristic energy exchanged 
in the electron-electron interactions scales with the in- 
verse of the dwell time Tl/td in each quantum dot - 
an observation that is closely related to the uniformity 
of the interaction potential inside a quantum dot. The 
number of quanta exchanged along a typical trajectory is 
too small to lead to a significant broadening of the distri- 
bution function — in that sense transport in a quantum 
dot network is always quasi-elastic — , although the ex- 
change of a single quantum is sufficient to suppress the 
interference from time-reversed trajectories. 

The semiclassical arguments of this section relied on 
the treatment of <j>(r, t) as a classical fluctuating poten- 
tial. In this respect, we follow earlier works on quantum 
dots by Seelig_and ButtikerEa and on granular metals by 
Blanter et aZ.Eil This approach was taken originally by 
Altshuler et al. for dephasing in quasi oae-dimensional 
and two-dimensional disordered metals. Q In the next 
section, we confirm the validity of this approach in the 



present context by performing a fully quantum mechani- 
cal calculation of the weak localization correction to first 
order in the interaction propagator L. The calculation 
of Sec. IV shows that the potential fluctuations are es- 
sentially classical if T > h/ro, where td is the (typical) 
dwell time in a quantum dot in the network. Since H/td 
is much smaller than the relevant temperature scale 
for the suppression of the weak localization correction by 
electron-electron interactions, cf. Eq. (^4|) of Sec. Q, this 
proves the validity of our approach for all temperatures 
of interest. 



IV. QUANTUM MECHANICAL CALCULATION 
A. Random matrix formulation 

We consider a network of Afo chaotic quantum dots 
coupled to electron reservoirs. The Hamiltonian of the 
entire system is written as 



H = H o + H u 



(43) 



where Hq describes the electrons inside the quantum dots 
or inside leads without taking into account their interac- 
tions, and £/i n t describes the electron-electron interac- 
tions. We write the non-interacting Hamiltonian Hq as 
a sum of three terms, 



Hi fl — Ht 



H 



DL 



(44) 



where Hd and Hi, describe the electrons inside the 
quantum dot network and inside the leads, respectively, 
whereas Hm describes the coupling between the quan- 
tum dots and the leads. We now describe each of the 
three terms contributing to H separately. 

Linearizing the electronic spectrum around the Fermi 
energy inside the leads, we have 



a=l,2 j=l 



dk 

2^ 



VajWlAk)il> a Ak), (45) 



where the index a = 1,2 labels leads connecting to 
the left and right electron reservoirs. The operators 
VajCO an d j(^) are for electrons in scattering states 
at wavenumber k (measured with respect to the Fermi 
wavenumber) and transverse mode j. The total number 
of propagating modes in the leads connecting to reservoir 
a is N a , a = 1,2. [If a reservoir is coupled to more than 
one lead, the summation over the index j represents a 
sum over the transverse modes in all leads connected to 
the given reservoir.] Finally, u Q j is the Fermi velocity of 
electrons in mode j. The current operator I a reads 

N a 
J'=l 

(46) 
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where 



dk 
— < 
2tt 



,±ik8 



1> a Ak), a =1,2, 



(47) 



where the iVj x iV a matrices Wj a = are related to 
the transmission matrices tj a of the contact between the 
ith quantum dot and reservoir a, 



and S > is a positive infinitesimal. 

We use random matrix theory to describe the quan- 
tum dots. Following standard procedures, the electron 
operators in each quantum dot are represented by an 
Mj-component vector ^ . , where the index j = 1, . . . , A/d 
labels the quantum dots in the network and Mj is the 
dimension of the subspace corresponding to the dot with 
index j. The Hamiltonian then reads 

AT D Mi 

tin = E E 4t«Bi,aAe 

i=l «,j9=l 



(48) 



i<3 q,/3 



Here the elements Hi^ a p of the Mi-dimensional matri- 
ces flj are random numbers taken from from a Gaussian 
distribution with zero mean and with variance 



(-Hi, a/3 -Hi, 7<5) 



(Hi.apHls-y) 



A' 



M, 



SaSdpy + -^jrSaySpS- (49) 



t ia = 2nW ia (u a u i M i ) 1 /\M i + n 2 u i ^ 2 W ai W ia u 1 J 2 )- 1 , 

with a = 1,2 and v a is an 7V a -dimensional matrix with el- 
ements {v a )ij = Sij t(2irhv a> j) . The dimensionless con- 
ductance Qij and and form factor /y of the contact be- 
tween dots i and j are defined in terms of the transmis- 
sion matrix ty as in Eq. (^J) . Similarly, the dimensionless 
conductance g' ia — g' ai and form factor f' ia — f' ai between 
the dots and the two electron reservoirs are defined in 
terms of t' ia as in Eq. (bh. 

For the electron-electron interaction we take density 
fluctuations inside each dot to be well screened, so that 
the interaction couples to the total charges of the dots 
only, 



E 



-m 



Mi 

E 

a=l 



The parameters Xi and A^ are related to the density pf 
states V{ and magnetic flux <j>i in each quantum dot, ell 
i = l,... ,Mr>, 



Xi = 



M 



X' = 



M 



1 - 



4A/. i $2 



(50) 



where $o the flux quantum and £rh,i is the Thouless 
energy of the ith quantum dot. Further, in Eq. (48) 



the Mi x Mj matrices Vij are related to the transmission 
matrices ty of the contact between dots i and j, 

Uj = 2nV ij {viVjMiM j ) 1/2 (MiM j +TT S v i v j vl j V ij )- 1 . 

The Hamiltonian ifoL describing the coupling between 
the dots and the leads reads 



2 N a Wd Mj 

a* = EEEE 

a—1 j — 1 i—1 a—1 



dk 

27 



(52) 



where the capacitance matrix C was defined in Eq. ( |l2| ) 
above. The corresponding interaction Hamiltonian for 
a single quantum dot is known as 'universal interaction 
Hamiltonian'. c3 

Evaluating the conductance g of the quantum dot net- 
work and its leading interaction corrections using the 
Kubo formula one finds 



G = 



-g, g = go + 5g de v h + 



(55) 



where go is the conductance in the absence of interac- 
tions (i.e., for Hamiltonian Ho), and 8g dcph and Sg mt are 
interaction corrections. (The reason for the separation 
between 8g deph and Sg lnt is that these two corrections 
have different temperature dependences, as will become 
apparent later.) Denoting with "•" adjacent indices to be 
summed over [as in Eq. (|i"4|)], the three terms in Eq. (|5T 
read 



,9o 



in 2 J de [-dJ(e)}tYViW 1 .G R (e)W.2V2W 2 .G A (E)W.i, 



(56) 
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and the interaction corrections 5g deph and 5g lnt are 

Sg d ^ h = An 2 J de J ^ [-d s f(e)} [coth(u;/2T) + tanh((s - w)/2T)] £ Im [L*( w )] 



x tr 



iyiWi.G R l (e)Gf j (e ~ u)Gf.(e)W. 2 v2W 2 .G A {e)W. i 



+ v 1 W 1 .G R (e)W. 2 v 2 W 2 .G A i (e)G A {e - u;)Gf.(e)W. 1 



+ -ujWt.G^e ~ u)Gf:{e)W. 2 iy 2 W 2 .G A {e)Gf.{e - w)W 



+ - l / 1 W 1 .G^(e)G?'(£ - uo)W. 2 u 2 W 2 .G A \e - u>)G A {e)W. 1 



(57) 



hit 



4^ 



2^ 



.Vd 



-9 £ /( £ )]tanh((e-c)/2T) ]T Im|L&(w) 



x tr [i/ 1 Wi.G^(e)G?-(e - w)G*(£)IT. 2 ^ 2 IF 2 .G A (£)VFi 

+ ^i^i.G R (£)^. 2 ^ 2 ^ 2 .G^(£)G^(£ - w)G£(e)Wi] 



(58) 



In these equations G^ and G A ; denote the retarded and 
advanced Green functions of the network of quantum 
dots without the electron-electron interaction Hamilto- 
nian H\„ t . These are matrices of dimension M, x M. 



lint 

which are the solution of 

2 

£ — Hi + in 

a=l 
2 

£ - Hi - ITT ^ WiaVaWa 



J' 



a=l 



Gg( e )+V i .G3( E ) = ti 

G A (£) + K.G A (£) = 1, 



(59) 



(H af ,H^ s ) = 



A yrf bb + A rfrf b'o 

a/3 7 5 a/3 




with lj the Mj x Mj unit matrix. Finally, L^(w) and 
L A (w) = L^(w)* represent the (RPA) screened interac- 
tion propagator, see Eq. (|3|) above. 

It remains to calculate the ensemble average of the con- 
ductance G for the ensemble of Hamiltonians described 
by Eq. ([l9|) above. This is the subject of the next sub- 
section. 



FIG. 4: (a) Diagrammatic rules for the ensemble average us- 
ing Random Matrix Theory. The weight factors depend on 
the symmetry present: A' = A in the presence of time rever- 
sal symmetry, while A' is reduced in the presence of a weak 
magnetic field and A' = where time reversal symmetry is 
fully broken, (b) Expansion of the full matrix propagator in 
terms of single propagators \/{e-\-mvWW^), depicted by sin- 
gle lines, and the matrix elements H a/ s, depicted by two open 
circles, (c) Dyson equation for the self energy E. 



B. Average over random Hamiltonian 

The average over the random matrices Hi is per- 
formed using a .-variation of the impurity diagram- 
matic technique JIfl This technique has been applied 
for various transport and thermodynamic properties 
of chaotic . quantum dots without electron-electron 
interactions.E3E3E3L3 Below we present its generalization 
to arbitrary networks. 



1. Average Green function 

We first discuss the calculation of the ensemble average 
of the Green function, (G R (e)) and (G^fe)). Following 
the diagrammatic rules laid out in Fig. |4j and keeping 
diagrams in the non-crossing approximation only,L3 i.e. 
diagrams without crossing double lines, one finds that 
the ensemble averaged Green function (Gfj(e)) satisfies 



11 



b) 



— «-Cp <p — — Hp (p— *(p Ij>— *■ — "-Ip lp— *(p tp— *(p <p* — 



— *-c> rj-* — *q^t^-*p p-*- — »»qc}j*c> 9— o— * 

= •. + + ;>fC v h 

»-<p p — ► — *<p p~~*p <p — * — Hp p~^p <p~~*p (p — * 

= : ; + MM + MMM -I 

Jx-»-6 ^a y* b-^- 1 ^ O > 6 J)s , 



FIG. 5: Diffuson ladder (a) and Cooperon ladder (b). 

the Dyson equation 

(G§(e)> - G?(e)« +X)G?(e) ifc E fc <G^(e)>, (60) 



where the self energy is 



E*( £ ) = -ftr(G*( £ )>, 



(61) 



and G^(e) is the solution of Eq. fl59| ) with rT^ = 0. Com- 
bining Eqs. ( |60|) and (|6l]) gives a self-consistent equation 
for E R . In the limit M{ 3> <7- x + g^ 2 + 5y > one finds 

(G?-( £ )) = (G^( £ ))t 



J7T 



Mi + Aj 

7T 



AMiMj 



2Mj 



■kv;E — ztr 



A,- 



M» + Aj 



%,(62) 



where i'y and ijj are given in Eq. ( |11| ) and (|5lj) above 
and A, is an hermitian M; x M; matrix, 



Nn j 2 

A, = tt 2 ^ V— Vife^feT4 4 + 7r 2 ^ yV ia z; a W ai . (63) 



a=l 



2. Classical conductance 

To leading order in the average number N of trans- 
mitting channels per dot, the calculation of the average 
conductance involves the calculation of geometric series 
involving the ensemble averaged Green functions. Dia- 
grammatically, these geometric series correspond to "lad- 
der diagrams" , as shown in Fig. || Such ladders are the 
equivalent of the "diffuson" propagator in diagrammatic 
perturbation theory. The building block of the geometric 
series is 



tr(G*(e))<G£(e')> = 



Mi 6lJ 



(64) 



4MM, 



[g - i2ir(e - e')v\ 



a) 



b) 



> b <MM-y 

A 

d) 




FIG. 6: (a) Diagrammatic representation of the leading con- 
tribution g cl to the ensemble-averaged conductance (g). (b)- 
(e) Diagrams contributing to the weak localization correction 
5g wh . (f) Definition of the Hikami-box used in (c)-(e). 



where g^ was defined in Eq. (Q) above. Summing the 
geometric series in Fig. ||a then gives the Diffuson matrix 



Dij(e,e') 



2Mi 



[g — i2tt(e — e')i> 



V 7T1 / ■ ' 



(65) 



For the calculation of the mean conductance one also 
needs a trace that involves the lead indices, 



D' ia = nv a tr [W ai {Gl){G±)W ic 



4M 



a = 1,2. 



(66) 



Combining everything, we then find the leading conduc- 
tance of the system 

(g)=g[.(g- 1 )-9'.2, 

which is equation (|lj) of section 0. 

3. Weak localization correction 

The above calculation gives the conductance to lead- 
ing order in g. A correction to sub- leading order in g is 
given by a class of diagrams that contains a maximally 
crossed ladder, as shown in Fig. ||b. These contributions 
are analogous to the "Cooperoa" contributions in dia- 
grammatic perturbation theory.EJ The summation of the 
geometric series promotes the contribution to be of or- 
der 1/N instead of 1/M, as is the naive expectations for 
diagrams that contain one crossed line. 
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5{G) 



J6= 




As this contribution is already small as 1/M, one may 
neglect the effect of a weak magnetic field on this term. 
The self energy correction <5£ affects the diffusion ladders 
as D -> D + 5D, with 



FIG. 7: Dyson equation for corrections to (Gu) due to the 
possibility of Cooperon like ladders in the time reversal sym- 
metric case. Double-hatching indicates a retarded-retarded or 
advanced- advanced pairing. These ladders are parametrically 
small, and for that reason can also not extend across multiple 
dots. 



2M. 



2 U V 



(70) 



In contrast to the diffuson propagator discussed above, 
the cooperon propagator is sensitive to magnetic flux. 
Proceeding as before, we hnd 



Cjj (e, e ) 



-[ga+g- 2wii>(e - e')}. 1 — ^ (67) 



with defined in Eq. ([l0|). For the calculations below, 
we also need geometric series of Green functions of the 
same type. These read 

Cf^e') = Of ( e , £ ')* 

= 77r4 f ( 8M i + ffH,« + <7ii 

- i2ir(s + e>i)*« - hi (1 - %) ) • (68) 

Cooperon ladders give a correction to the self-energy 
appearing in the calculation of the average Green func- 
tion, as depicted in Fig. ^. Calculation of the self-energy 
correction SY,i to leading order in g/M then gives 

A,- 



= ^tr[(G*)(Cr<G*) + ^)<G*} 



(69) 



This contribution is depicted in figure ^|b. 

In the diagrams for the weak localization correction 
to the conductance, the Cooperon and diffuson propa- 
gators are connected in a so-called "Hikami box" JJJ In 
our diagrammatic analysis the analogue of a Hikami box 
is depicted in figure ||f. We consider the general case of 
a Hikami box with four energy arguments. We write E\ 
(e[) for the energy argument of the retarded (advanced) 
matrix propagator on the left side, and £2 (e' 2 ) for the 
energy argument of the retarded (advanced) propagator 
on the right. For the calculation of the weak localiza- 
tion correction one only needs the case of equal argu- 
ments, £\ = Bi = £2 = s' 2 - For dephasing and interaction 
corrections, some arguments differ. Explicit calculation 
shows that the Hikami box depends on the combination 
u) = e[—ei+e' 2 — £2 only. Hence we write Bij t ki(uj), where 
the indices i and j refer to the left and right (Diffuson) 
ladders and the indices k and I refer to the bottom and 
top (Cooperon) ladders. 

The calculation is essential but technical; we outline 
it in the appendix . The Hikami box Bij t ki(to) is zero 
except where at most two different indices appear, 



Bij : ki (w) 



IQMiMiMkM, 



2iriviUj8ij5j k 5M + (5uS jk + SikSji)(gij + m,ij - fij) 



+ (dikSu + Sjk$ji)fij + {SijSki + SijSu)fki — SijSkifik 



(71) 



For the evaluation of the weak localization correction, one also needs to consider Hikami boxes that are connected 
to the leads, not only to Diffuson propagators inside the quantum dot network. The two contributions of this type 
are depicted in figure ||c and d. They are 



B' 



7T 3 Z/? , 
— f . 

16Mf Jar 



(72) 
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FIG. 8: Diagrams for the first-order dephasing correction. Diagrams depicted in (b), (c) and (e), (f) are weighed with a factor 
1/2, in line with Eq. (|57|). Together (a), (b) and (c) constitute the correction to the Diffuson propagator, which cancels to 
leading order. Hence the only relevant contributions are the corrections to the Cooperon in (d), (e) and (f). In both cases, 
complex conjugate contributions exist which are obtained by placing the vertices on the opposite matrix propagation lines. 



Combining everything, we have (see Fig. ra) 



No 

Sg WL = 2D' 1 D..8D..D..D[ 2 • 2 ^ ( ", , (D[D..B'. 2 jl + lt J)A)' 2 + D' 1 .D..B.. tji (0)D..D'. 2 ) , (73) 

where D' ia = D' ai was defined in Eq. ( |66| ) above and we have suppressed superscripts as well as inconsequential energy 
arguments of Z? RA (e,e), C RA (e,e), cf. Eqs. (|65|), (|68|). The four terms correspond to the four diagrams b - e of Fig.||. 
Substituting our results for the Hikami box B, the Cooperon and Diffuson propagators C and D, and the interaction 
propagator L, we arrive at Eq. ( |l5| ) of Sec. y, with the zero-temperature Cooperon c = (g + <?h) _1 - 

So far we have not taken into account electron-electron interactions. To lowest order in perturbation theory in the 
interaction Hamiltonian -Hint, the dominant interaction correction to weak localization comes from Sg deph in Eq. (|57|). 
The corresponding diagrams are depicted in Fig. |^. We now calculate that correction. This interaction correction 
is nonzero only if both interaction vertices appear inside the Cooperon propagator. (This is why this interaction 
correction does not affect the leading contribution go to the conductance.) 

To calculate the interaction correction, one notices that the interaction vertices are "dressed" , as is shown in Fig. 
HI For this case energy arguments may be neglected, as they lead to corrections small in g/M. Labeling the dot in 
which the interaction takes place by the index a, the dressed interaction then reads 

7 R .. = (l A ..)* 

= tr [<G R ) (1 + tr [(G R )(G R )] D™) (G R )(G A )] 6 ai 6 aj 

2Mi 1 31 2Mi V ' 



where 



1 r -| 

(8Mi + g u - i2n(e + e')^)% ~ 9tj C 1 ~ %) • (75) 



The interaction correction SC to the equal-energy Cooperon propagator G(e,e) then becomes 



SCa = 



JdsJ ^[-S e /( E )][coth( W /2T) + tanh((e- W )/2T)] g Tm[L%(u)] 



a,/3=l 



x [G.(£,£)/ R ..G..(e-w,e)/ R ..G. J (£,e) + G.(e-w,£)/ A : ..G..(e-w,£-w)/ R ..G. J (£,£-^ . (76) 
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FIG. 9: Renormalization of the interaction vertex by ladder diagrams involving Green's functions of the same type (retarded- 
retarded or advanced-advanced). 
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FIG. 10: Dyson equation for the Cooperon obtained by perturbation theory in the high temperature limit. The hatched boxes 
indicate noninteracting Cooperon ladders, while gray shading indicates that interactions are taken into account. Wiggly lines 
indicate the equal time interaction propagator, which can either connect back to the same propagation line, or to the opposite, 
time reversed one. 



Performing the energy integration and passing to dimensionless propagators, we then find 

Sea = f ^ % V* Im \^ 2 v a vnL^Juj)\ (77) 

3 J 2ir 2T sinh 2 (w / '2T) ^/^v V ; 

x {(<? + <?h + i27rwi>)7 Q 1 (g + g H )-* (g + g H - i2iruD)^j - (g + g H )^ (.9 + 9n + i27rwf?)~^ (g + gn)^] 

+ (.<? + m - i2nuji>)^ (g + H )-* (g + g H + i2'xuv)~ 1 j -{g + gn)~a (g + m - i2nwv)~p (g + } - 

Let us now inspect the integral in Eq. (|77|). The term between brackets {. . .} is proportional to lj~ 2 if uj > ?i/td, 
where Tl/tv ~ g/v is the inverse dwell time of a dot in the network. Since ImL R (w) oc u> for u> ~ h/m, one thus 
concludes that the integral in Eq. ( ff^ ) converges at uj ~ min(?l/TD, T). We focus on the regime T ^> Ti/td, in which 
the convergence is at cj ~ U/t^. In this regime the inequality w < T is obeyed for all frequencies w contributing to 
the integral, so that all relevant interaction modes that contribute to dephasing can be described using the classical 
fluctuation-dissipation theorem. Indeed, one verifies that in this regime the first-order interaction correction (]77 



agrees with the interaction correction to c obtained in the semiclassical framework of Sec. Ill, taken to first order in 
the interaction propagator L. 

Estimating the magnitude of the first-order correction Scij for T 3> Ti/td 1 we find that Scij ~ CijT/T^, where 
~ hg/TD [see Eq. ( p4| ) above]. This observation has two consequences: First, it implies that the regimes of 



validity of first-order perturbation theory and the semiclassical approach of Sec. Ill overlap: Both approaches are 
valid if H/td <C T <C T^. Second, it implies that interactions give no significant correction to the weak localization 
correction 5g WL if T < SVtd, so that we may ignore the difference between the fully quantum-mechanical interaction 
correction S£ij of Eq. ( yT\) and the semiclassical result in the low-temperature regime T < fi/r^ within the limiting 
procedure outlined in Sec. H (Both approaches give essentially no interaction correction to weak localization at these 



temperatures.) When combined, these two observations justify the semiclassical considerations of Sec. Ill, as well as 
the expressions ( |l5| ) - ([l?]) for the weak localization correction Sg wh that followed from these considerations. 

For completeness, we mention that the full temperature dependence of 5g WL can also be obtained from diagrammatic 
perturbation theory. Following the above arguments, in the limit T ~^>%/td all factors coth(o;/2T) + tanh(e — ui)/2T 
appearing in the calculation may be replaced by 2T/u>, irrespective of the value of e. This considerably simplifies the 
calculation, and the m interaction propagators that appear in mth order in perturbation theory may then be placed 
independently of each other along the cooperon ladder. Using Eq. (p9|) for the interaction propagator and writing the 
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Cooperon ladders (without interaction corrections) in an integral form similar to Eq. ( J3S 

/>oo 

(g + gn + ^muv)" 1 = (2ttTH>)~ 1 / dTe-^~ l " r , (78) 



one may perform the frequency integrations. The resulting expression consists solely of time integrations with instan- 
taneous interactions. The remaining combinatorial problem leads to a Dyson equation of the form shown in Fig. |Tc| . 
Here the first term on the right hand side is the noninteracting Cooperon c^i = (g + 9Yi)ki and the six other terms are 
obtained by different placements of the interaction propagators. [Note that where beginning and end are on the same 
Green's function line, an additional weight of 1/2 arises from a factor f?° cIt8(t) = 1/2.] Adding the six different 
contributions gives a vertex proportional to {AitT / d s fi)(g~^ m +g^n ~ 2 9mn)^ so that one arrives at the Dyson equation 



Mo 

<id = (g + m)ki 1 - Y] (r + r H ) _1 r s mn , (79) 

— _ L J kmAn 



m,n— 1 



where T, Th, and are rank- four tensors whose definition is given below Eq. (0). With a little algebra, one verifies 
that Eq. ([79|) is equivalent to the result (1(J) derived using semiclassical arguments. 

Equation ( fzr| ) can also be used to calculate the magnitude of energy quanta uj exchanged with the fluctuating 
electromagnetic field in the quantum dots. Hereto, we note that the sum of the second and fourth terms between 
brackets {. . .} in Eq. ( |77| ) is proportional to (minus) the probability pi(u>) for emission or absorption of a photon 
along the electron's trajectory, so that 



PlM = , — i , o rr ■ i2 ^ M^V^a/jMlMfll-ff-a 1 (fl + i27rUji) )a0 9^ 9-2 

9i-9- 9-2 £pL x 27rTsmh w/2T 



16Tw 2 



7V D 

Y u «9ap v P^Wi-9^(§ + i^i>)~pg^g' 2 ], (80) 

a, (3=1 

where, in the second equality, we took the limit T ^> h/rp. The probability that one inelastic scattering event of 
arbitrary frequency occurs is Pi = j dujpi{uj). Equation ( |80| ) is valid as long as Pi <C 1, so that first-order perturbation 
theory is sufficient. 

From Eq. ( |8C| ) we conclude that the energy of photons that are emitted or absorbed is limited by mm(h / td , T) . 
The temperature at which the interaction correction to weak localization becomes relevant is the temperature at 
which the probability that at least one energy quantum is exchanged becomes of order unity. However, the typical 
exchanged energy remains of order h/m for all temperatures. This implies that the broadening of the distribution 
function by inelastic processes is parametrically smaller than the temperature T, by a factor 1/g <C 1. Transport 
in the quantum dot network is thus quasielastic for all temperatures. (Inelastic processes become relevant only if 
T > pThjiff 1 / 2 , where -Erh.z is the Thouless energy of an individual quantum dot.) 



4- Interaction corrections to the conductance 



The relevant diagrams for the interaction correction to the conductance <5g mt are shown in Fig. [H]. These diagrams 
do not involve Cooperon propagators. The diagram shown in Fig. [Ha is analogous to the ones we have already 
encountered in calculating the (first-order) dephasing correction to weak ocalization. It gives an interaction correction 
to the diffuson propagator D(e,e) that depends on the frequency to of the interaction propagator, 

5Dp atij (u)M = D i .(e,e)lf t ..D..(e-u,e)I*..D. j 

= " — 9ip (9 + iZkuv)^ g aj — ■ (81) 

(The frequency to will be integrated over in the final expression.) For the remaining diagrams, we need to consider an 
interaction vertex that connects an advanced and a retarded Green function. Such an interaction vertex is dressed by 
a Diffuson propagator, which allows the interaction vertex to be placed in a dot different from the one that appears 
at the outer end of the dressed interaction vertex, 

= S m + Y / D lk (e-u J ,e)ti-(Gt(e))(G^ k (e-u J )) 

k 

4M^ 



V; 



(g + i^ujv)^ . (82) 
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FIG. 11: Diagrams contributing to 5g mt . The Hikami box is denned in Fig. ^[ 



With this interaction vertex, the diagrams of Fig. [Tlp- d (without the outer diffusion ladders) can be represented by 
Hikami boxes Bij^iiyj) and B' a ^ kl of Eqs. ( |7l| ) and (|72|), but with — > because no Cooperon ladders are involved. 
Combining the two interaction contributions to the interaction correction we find 



-,int 



dui f d 



2-7T V du> 



a/3 k,l=l 



L%(w)I k J w (B' llrk D..D'. 2 + D[D..B[ l>2k + D' 1 D..B. l , k (uj)D..D[ 2 



(83) 



Expressing the propagators in terms of the matrices g 
and /, we find that 6g lnt naturally separates into two 
contributions, which are given by equations (|20|)-(p2|) of 
Sec. ||. Both corrections are small for all temperatures, 
and it is not necessary to consider higher order contri- 
butions involving more than one interaction propagator 
L. 



V. APPLICATION TO DOUBLE QUANTUM 
DOT 

We now apply the theory of the previous sections to 
the case of a double quantum dot. There are two cases 
of interest: A linear configuration, in which each dot is 
coupled to one reservoir, see Fig. [l^a, and a side-coupled 
configuration, in which both reservoirs are connected to 
the same quantum dot, see Fig. [T^b. 



A. Linear configuration 

The conductance matrix for the linear double quantum 
dot reads 

9 = ( 9 ' 11+9U rT ), (84) 
a V -912 g' 22 + 912 J y ' 

where g' lx and g' 22 are the dimcnsionlcss conductances of 
the contacts connecting the two dots to the reservoirs, 
and g\2 is the dimensionless conductance of the contact 
between the two dots, see Fig. [l2[ The form factor matrix 
/ has a similar structure, with g[ ly g' 22l and g\2 replaced 



by /id /22> an< 4 /i2j respectively. The classical conduc- 
tance of the system is G c \ = (d s e 2 /h)g c \, with 

.9 c l 1 =.9n 1 +.9 2 2 1 +5r 2 1 , (85) 
see Sec. g, Eq. (||). 

1. Weak localization 

The zero temperature weak localization correction to 
the conductance 5G WL = (d s e 2 /h)Sg WL follows from 
substitution of the zero-temperature Cooperon c(0) of 
Eq. © into Eq. @, 

Sg WL = f[j9'A + h2/gl2 

9 2 c \ 9'u + 3h : i + gi2 - 5? 2 /(.9 2 2 + 5h,2 + 512) 

/22/ff22 + f 12 1 9 12 

g'22 + 9H,2 + 312 - 9l 2 / (9'u + 9H,1 + 312) 
2(ffl2 - Zl2)/.9l2 

(9'u + 3H4 + 9i2)(3 2 2 + 9H,2 + 312) - 3l2 ' 

(86) 

Here (7H.2 and gu.i are dimensionless numbers describing 
the effect of an applied magnetic field, see Eq. (|Io| ) . The 
limit of zero magnetic field c/h,2 = 3h,i = agrees with, 
the result obtained previously by Golubev and Zaikin.lEl 
The high-temperature limit of 5g WL ' d of the weak local- 
ization correction is found by taking the diagonal contri- 
bution c d of Eq. ( |l9| ) for the Cooperon propagator, 

5g WLA = f'n/g'A + fu/gh Wg^+^/gj 

3ci g'n + 3h,i + 312 g'22 + 3h,2 + 312 

(87) 
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FIG. 12: Schematic drawings of two double quantum dots. 
Panel (a) shows a linear configuration; Panel (b) shows a side- 
coupled configuration. 



Note that |<5<?wlI < |^<7wl|- The remainder of the weak 
localization correction, Sg WL —5g WL,d , is temperature de- 
pendent because of dephasing from electron-electron in- 
teractions. Taking the temperature-dependent Cooperon 
from Eq. ([jl]), we find that the temperature dependence 
of the full matrix c(T) is encoded in a single scalar func- 
tion/(T), 



c(T) = 2(0) -[c(0)-c d ]/(T). 

Equation (pq) immediately implies that 

5g WL (T) = 6g WL - d + [6g WL (0)-6g WL - d ][l~f(T)}, (89) 

where 5g WL (0) and Sg wh,d are given in Eqs. ( |S6| ) and 
(|S7|), respectively. In the regime where temperature is 
large enough for dephasing effects to give a sizeable cor- 
rection to the weak localization correction to the conduc- 
tance, we obtain f(T) from Eq. (|l6|), 



f(T) = 



T 



with 



T± = Hn + T 2 )(g' 11 g' 22 + g' n gi2 + g'22912) 
d s 4jrr + T-(g' n + g' 22 ) 



(90) 



(91) 



Here t\ and t 2 are the (classical) dwell times of the two 
dots, modified for the presence of a magnetic field, 



ti = 



9n + .9H,i + gi2 



T 2 = 



2TtHl>2 



9 22 + .9H,2 + 312 



-,(92) 



whereas t± are time scales representing the relaxation of 
symmetric (+) or antisymmetric (— ) charge configura- 
tions in the double dot, 



(93) 




n 2 h 2 v\v 2 

It is instructive to compare Eq. (KXJ) with the expres- 



sion for f(T) obtained in first-order perturbation theory, 
duj lo/2T 2uj 2 



f(T) 



2-k sinh 2 (w/2T) (1 + Lu 2 Tl/h 2 ){\ + uj 2 t 2 /% 2 ) 
n t\t 2 



Im[L? 1 (u)+L*(uj)-2LK(L 



(94) 



The integral in Eq. (Q) converges for frequencies lo/Tl 
of order t± . For these frequencies, we may neglect the 
capacitance C in the expression for the interaction prop- 
agator L since C/e 2 <C v. The resulting frequency inte- 
gration yields 



f(T) = 



2-kTt+t- T 
3H(t + - t-) T$ 

x \T x (2txTt-IK) - J- 1 (2ttTt +/h)] (95) 



where 



2ip' 



(96) 



tp' being the derivative of the digamma function. With 
the asymptotic behavior of J-\(x), 



1 2 

5 X 



X < 1, 
X > 1, 



(97) 



we identify three different regimes for the temperature 
dependence of the dephasing correction: 



f(T) = -^T + T_(r + 



. /2nT\ 3 T , s 

T -» hr % < 98 » 



if T < h/r+ 



f(T) 



2ttTt_ T 
3h Ta 



if H/t + < T < h/r-, and 



f(T) = T/T^ 



(99) 



(100) 



if Tl/t^ < T, where is given by Eq. @ above. The 
intermediate temperature regime exists only if r+ ^> t_ . 
A comparison of Eq. ( |100p with Eqs. ( p0| ) shows that 
the two expressions for f(T) agree in the temperature 
regime Ti/t- <C T <C T$ where both expressions are 
valid. It is in this temperature regime that the factor 
(w/2T)/sinh 2 (w/2T) in Eq. @ can be approximated 
by 2T/lo, which is the appropriate weight appearing in 
the classical fluctuation-dissipation theorem. 

It should be noted that the low temperature correc- 
tions, Eqs. ( |9g| ) and (|99|), result in contributions to the 
conductance of order 0(1/ g). Such contributions are 
beyond the accuracy achieved in the limiting procedure 
outlined in Sec. U Further contributions of the same or- 
der might be obtained by calculating, e.g., weak localiza- 
tion corrections to the interaction corrections Jg 111 *' 1 and 
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5g mt ' 2 . For disordered metals such contributions have 
been considered explicitly in Ref. fj]. 

The above equations take a simpler form in the limiting 
cases of large and small interdot coupling 312 and of a 
large magnetic field. For small interdot coupling 1712 <C 
mm(9ii = 922)1 one has 



WL 



/i25n /22512 + /12522 



ffnKi + 3h : i) 522(522 + 



2(312 - /l2)ff!2 



3ll322 



(flii + 9n.i)(g' 2 2 + m,2) T^ + T' 
{9'xi + dn,!)"! 1 + (.922 + 9n,%)v2 l 



(101) 
(102) 



so that only a small part of the total weak lo- 
calization correction is temperature dependent. In 
the limit of a large interdot conductance, 312 3> 
max(3^, g'22, 9h,i, 311,2), the full weak localization correc- 
tion acquires a temperature dependence, 



WL 



922/11 + 911/22 



(g'll + 322) 2 (3ll + 3H,1 + 922 + 9H,2) 



312 



T + T' 

{g'll + 3H,1 + g'22 + 3H,2)(^f 1 



8ir 2 



-(103) 



Finally, in the limit of large magnetic field, 3h,i,3h,2 3> 
max(3 11; 3 22 ,3i2), we have 



Sg 



WL 



2 fill 9ll + J "l2/9l2 2 Z22/922 + i 'l2/9l2 

9cl — 



3H,1 

2 912 - /l2 
~~ 9cl 



9H,2 



312 
8tt 2 



3i23h,i3h,2 T^ + T 
[gn.1^1 1 +gn,2V 2 ~ 1 ) 



(104) 
(105) 



A special case of two weakly coupled quantum dots 
(312 < g'n, g'22) witn tunneling contacts {f' xl < g' n , 
I22 ^ 322, /12 <^ 312^-ihas been considered recently 
by Golubev and ZaikinS While our calculation agrees 
with that of Ref. ^33] in the high temperature regime 
T T$, significant differences appear in the low tem- 
perature limit. In particular, Golubev and Zaikin find a 



finite dephasing correction to weak localization at zero 
temperature, whereas we find no such effect. A similar 
discrepancy has been found previously in the context of 
dephasing froj»-the electron-electron interaction in disor- 
dered metalsJ3e3 In this case the neglect of recoil effects 
in the influence functional approach used by Golubev and 
Zaikin has been identified as the cause of the problcm.Eil 
This causes an ultraviolet divergence, which does not ap- 
pear in the perturbation theory, where it is avoided by the 
tanh-term in the factor coth(w/2T) + tanh ((e — u>)/2T) 
that sets the magnitude of the dephasing correction at 
low temperatures, see, e.g., Eq. (|57]) and Refs. |7|,|Sl|. (Ne- 
glect of recoil amounts to neglecting the w-dependence of 
the argument of the tanh, which causes this factor to no 
longer approach zero at large frequencies cj.) We believe 
that the discrepancy between our result and that of Ref. 
|63| has the same origin. 



2. Interaction corrections 

The interaction corrections Sg 111 *"' 1 and 8g mt 2 do not 
depend on the magnetic field. Hence, the relevant time 
scales do not involve 3h,i and 3h,2, and we define 



g'n + 9i2 



i = 1,2. 



(106) 



Again, we introduce time scales t± related to n and T2 
as in Eq. (^) above. For the first interaction correction 
jgint,! we then find 



int.l 



a 3 
9 c i 



d s3n3i23 2 2 



d uj 
—- UJ coth — — 
dw 2T 



xlm 



(r + +r_)/a 



(1 + iujT+/h){\ + itJT-/h) 



(107) 



This result was obtained previously in Ref. |3l] for the 
symmetric case g'n — g 22 , v\ = v% and in Ref. ^0| for 
the case g' xx = g' 2 2 — 312, vi — v-i- The frequency inte- 



gral in Eq. (107) can be evaluated in terms of digamma 
functions. We have 



J 



dio ( — — uj coth — I Im 
du 2T 



(1 + iujT a /Ti){\ + iujTpjK) 



2k 



T-2 



2irTT a 



-ft 



( h 



\2wTt0 



(108) 



where 



F 2 (x) = ip{l + x)+x^'{l+x) 



(109) 
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and 4>{x) is the digamma function.0 From the asymptotic behavior of T2, 



-7 + - 3C(3)x 2 + 



1 + lnx 



1 

12.r- 



X < 1, 
X > 1, 



(110) 



with 7 the Euler-Mascheroni constant, we obtain the high and low temperature limit of the interaction correction 



23c 3 i 



m — , r<7i/r±, 

T+ — T_ T_ 



d s9'ligi29' 2 2 



-kTi(t + + T_) 

6Tt + t_ 



T> ?i/r ± . 



(Ill) 



The second interaction correction <5<7 int ' 2 is expressed in terms of interaction-induced shifts Sg'^, 5g' 22 , and 8g\2 to the 
conductances </ n , g' 2 2, and 312, respectively, see Eq. (pl|). In contrast to the interaction correction ^g 111 *' 1 considered 
above, the frequency integrations needed to calculate Sg'n, 8g' 2 2, and 5g\2 converge only if we account for the finite 
(nonzero) capacitances of the quantum dots, see Eq. ( |22| ) . [The integration in Eq. ( |22| ) diverges logarithmically if the 
limit Cii/e 2 vi — > is taken.] 

Below we give explicit expressions for the case of a symmetric double dot only, g' xl = g' 22 = </ , f[% = J 22 = /') 
V\ = U2, and C — C\\ — 6*22- In this case, the logarithmic divergence of the integration in Eq. ( p2[ ) is cut off at the 
inverse of the charge-relaxation times, 



'c+ 



4eVC" d s e 2 v/{C + 2C 12 )' 

and the corrections £g' 1:L = <5g 2 2 — ^5' and are found to be 

f h \ „ / 7i 



55' 



^512 = 



.9' - f v 

2(512 



cr=± 



^2 



V27rTr 



/12) t + - r_ 



d s gi2 



T 2 



-T 2 
1 



2kTt cct 



2vrTr c _/?i 



(112) 



(113) 



(114) 



For the case (/ = (712, /' = /12 and C12 = 0, Eq s. (113) and (114) agree with results obtained previously in Ref. 
[3(3. [ The result of Ref. |3(1 differs from Eqs. (113) and (114) if C12 > because Ref. p0| includes cross capacitances 



between e ach d ot and adjacent reservoir of the same magnitude as the cross capacitance C12 between the two dots.] 
Equation |jj L 1 3 h s implifies to the renormalization of the contact conductance for a single quantum dot in the limit 
gi2 — * cxd.e3e3c3 Again making use of the asymptotic behavior of the digamma function, we find that the above 
expressions simplify to 



Sg' 



g'-f 

d s g' 



In 



x < 



In- 



,1+7 



Tr- 



,1+7 



2ttTt c+ 



I 6Tr + 



^ln- . 

r + 27rTr r _ 

T_ 



h/ T ± «T« h/r c ±, 
K/t c ± « T, 



(115) 



4(ffl2 - /l2) T+ - T_ 



d s gi2 



hi- 



x < 



7"+ 



In- 



e l+7 

27tTt c _ ' 



I 6Tr c 



t< h/r±, 

h/r- <T<C?i/r c _, 



(116) 



B. Side-coupled quantum dot 

For the side-coupled double dot configuration of figure |l2| the structure of the weak localization correction and the 
interaction corrections is essentially the same as for the linear configurations. The classical conductance is 

^i 1 =9n 1 +9i2 1 - (117) 
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The weak localization correction to the conductance is 
wl /22S11 + fn9i2 



Sg 



Ki + .9i 2 ) 2 (sii +9'i2 +9i2 +gn,i) i (g'n + g'12 + 9n,i)(9i2 + 9^,2) + gugn 



r 



where f{T) =T/(T <t> + T), 



1 Ti + T 2 

47T T+T_ 



-512, 



and 



2nhi>i 



g\i + 912 + m,i + gw 



To 



.912 + SH\2 



(118) 



-, (119) 



with r± given in terms of n and t 2 as in Eq. (|9£ 

Again, it is instructive to compare to what one finds to 
lowest order in perturbation theory. The result is iden- 
tical to Eq. (|95|), where ti , t 2 and T^> a re those of the 
side-coupled system, Eqs. (118) and (119). Simplified ex- 
pressions for the function f(T) in the regimes T <C h/T+, 
%/t+ < T < V 7 "-) an d ^/ T - < r are as in Eqs. (§|)- 



(100) 



In the limit of small interdot coupling g\ 2 — > only a 
very small fraction of the weak localization correction is 
temperature dependent, 



Sg 



WL 



/22.9n + /nSl2 



(.9ii +5i 2 ) 2 (5n +9i2 + .9H,i) 
1 



5 2 2 



Oil + .9i 2 + .9h,i).9h,2 T + T 



ff!2 

8tt 2 



g'12 + .9h,iK 1 + 3H\2^ 2 ■ (120) 



In the opposite limit of a large interdot conductance the 
entire weak localization correction is temperature depen- 
dent. In this limit there is no difference between the lin- 
ear and side-coupled c onfig urations, and one finds that 
Sg WL is given by Eq. (103) above, with g' 22 replaced by 
g' 12 ■ Finally, in the limit of large magnetic fields we find 



Sgwh 
d K 



22 



(9'n 



gi2 , 

8^ ^ 



g\2) 5h,i 
1 



1 



gh 



T lh 



gn.2V 2 



(121) 



With a side coupled quantum dot, the interaction cor- 
rection Jg 1 "*' 1 to the conductance vanishes. The interac- 
tion correction <5g lnt ' 2 coming from the renormalization of 
the contact conductances remains. The detailed expres- 
sions are rather lengthy and will not be reported here. 



VI. CONCLUSION 

We have calculated the quantum corrections to the 
conductance of a network of quantum dots, including 



WL 



the full dependence on temperature and magnetic field. 
Our results are valid in the limit that the quantum dot 
network has conductance g much larger than the con- 
ductance quantum, so that the quantum corrections are 
small in comparison to the classical conductance, and in 
the limit that the electron dynamics inside each quan- 
tum dot is ergodic. Following the literature, we sepa- 
rated the quantum corrections into the weak localization 
correction Sg WL and two interaction corrections Sg 111 *' 1 , 
Sg lnt ' 2 . Our results for the interaction corrections agree 
with previous calculations of Sg 11 ^' 1 and 8g mt 2 by Gol- 
ubev and Zaikinl23 for a linear array of quantum dots, and 
are closely related to similar interaction corrections in a 
granular metal, see Ref. E9l Our result for 5g 



with the literature in the limit o£-gero temperatureEro 
and in the high temperature limit ,Lil but we are not aware 
of a calculation of the full temperature dependence of 
5j WL in the literature. (The exception is a calculation of 
Sg WL for a double quantum dot by Golubev and Zaikin 
which, however, gpes an unphysical result in the limit of 
zero temperature.E3) 

We have formulated our final results in such a way that 
the evaluation of quantum corrections for a network of 
a relatively small number Ad of quantum dots does not 
require more than the inversion of an AD-dimensional 
matrix. All quantum corrections to the conductance can 
be expressed in terms of the inter-dot conductances, form 
factors, and the capacitances only. In principle, these pa- 
rameters can be measured independently. This makes a 
small quantum dot network an ideal model system to 
compare theory and experiment, and to unambiguously 
identify the mechanisms responsible for dephasing. (Ca- 
pacitances and form factors play a role only if the dots 
are connected via non-ideal contacts in which one or 
more transmission eigenvalues are smaller than one. For 
lateral quantum dot networks defined in semiconductor 
heterostructures, contacts are typically ballistic, and the 
only relevant parameters are the quantized conductances 
of the contacts between the quantum dots.) 

The simplest example of a small quantum dot network 
is a 'double quantum dot', which consists of two quan- 
tum dots coupled to each other and to electron reser- 
voirs via point contacts. Several groups have.re.port.erl. 
transport measuremepts on such double dots,EaHI£j - Efl 
or even on triple dots.E3 (Double quantum dots also play 
a prominent cole in recent attempts to achieve quantum 
computation.EJ However, the dots used in these experi- 
ments typically hold one or two electrons each and can 
not be described by random matrix theory.) Although, in 
principle, the contact conductances in lateral double and 
triple quantum dot networks are fully tunable, the exper- 
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iments of Refs. f2j||2||^,|84] were performed for the case 
that the devices are weakly coupled to the source and 
drain reservoirs. In that limit, transport is dominated by 
the Coulomb blockade. Our theory applies to the oppo- 
site regime in which all dots in the network are open, well 
coupled to source and/or drain reservoirs. We hope that 
the availability of quantitative theoretical predictions will 
lead to renewed experimental interest in quantum trans- 
port through open quantum dots. 
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APPENDIX: HIKAMI BOX CALCULATION 

In this appendix we provide details on the derivation 
of Eqs. (71) and ([72) of Sec. IV. The explicit expression 



for the Hikami box is an essential part of the calculation 
of the quantum corrections to the conductance, but we 
have not found the explicit expression of Eq. ( f7l| ) , nor its 
derivation, in the literature. 

We refer to the text surrounding Eq. (|7l]) for the nota- 
tions used in this appendix. In general, the Hikami box 
Bij,kl(u) will be nonzero only if the four indices span 
at most two adjacent quantum dots. We here show the 
calculation of Bh_u(lu). There arc three contributions to 
Bii,u(oj), which are shown in Fig. |l3[i,ii a-c. They read 



B {a) -(u>) 



tr [(GS( £l ))(G^(4))(G?|( e2 ))(G^( e ' 1 ))] 
mvi{s\ - e'i + £2 - 4) 



Mf 



2M,; 



tr 



-2Aj(Mj - Aj) Af 
(Mi + A,) 3 Mi(Mi + Aj) 4 
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2M 3 
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8Mi Mi 
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3MM, 
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= ( te [<G£( £ i))(G£(4))(G£( e i)]) G£ A (4,4) ( tr [(G£(4)><G£(e 2 )>(G£(4)]) 



2M 3 
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■ z2tt(2£i - 3e[ + 2e 2 - 3e 2 ) 1 



8M, 



Mi 



tr 



A 3 - 3MfAj 
(Mi + A,) 3 



(A.l) 



(A.2) 



(A.3) 



where the Mi x Mi matrix Aj was defined in Eq. ( |63] ) 
above and u) = e'-y — £\ + e' 2 — Traces involving the 



matrices Aj can be calculated using the identities 



tr 



tr 
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FIG. 14: Diagrammatic depiction of contribution from 
Hikami boxes placed adjacent to leads. 



for i 7^ j. Other contributions are related by symmetry. 
Rewriting the general case Bij t ki(u>) in terms of the ma- 
trices g and / for contact conductances and form factors, 
we obtain the result given in Eq. (|7l|) of Sec. IV. 



If a Hikami box is placed adjacent to a lead, one finds 
the three contributions shown in Fig. [Ti]. Adding these 
we find, with the help of Eq. (A.5), 



Addition of Eqs. (A.l)-rtAJ) gives 



^4, ,4 r 

7T V, 



16M? 



2mviU> + 2(<?h,m + 9a) + fit (A. 6) 



The diagrams for the relevant contributions to 
Bij,ki(w) in which the indices differ are shown in the 
other panels of Fig. 13. Expressing these contributions in 
terms of the mat rices Aj an d pe rforming the traces with 
the help of Eqs. (|A.4|) and (A..5), we find 



Bij^j (w) 



Bu y ij(iv) 



WMfMf 



{.fij 9ij) j 



lQMfMj 

1 J 



(-fa) ■ 



(A.7) 
(A.8) 
(A.9) 



B' 



tr 
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(A.10) 



This is the result reported in Eq. (72) of the main text 
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